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Abstract 



We consider a free boundary problem for the p-Laplace operator 
which is related to the so-called Bernoulli free boundary problem. In 
this formulation, the classical boundary gradient condition is replaced 
by a condition on the distance between two different level surfaces of 
the solution. For suitable scalings our model converges to the classical 
Bernoulli problem; one of the advantages in this new formulation lies 
in the simplicity of the arguments, since one does not need to consider 
the boundary gradient. 

We shall study this problem in convex and other regimes, and es- 
tablish existence and qualitative theory. 

1 Introduction 

For a bounded domain K C M. N (N > 2), the well known (exterior) Bernoulli 
free boundary problem is to find a domain D K, with the requirement 
that the harmonic function u in the region ft\K with given boundary values 
satisfies a prescribed Neumann condition on <9f2; the boundary <9f2 is thus 
called Bernoulli free boundary. There is a vast literature on the subject, in- 
cluding by-now classical references [H [2j [U [131 EES • O ne may also replace 
the Laplacian operator by the p-Laplacian operator; this problem arises in 
various nonlinear flow laws, and several physical situations, e.g. electro- 
chemical machining and potential flow in fluid mechanics (see for instance 
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In this paper we propose a new formulation of the p-Laplacian Bernoulli 
problem for 1 < p < +00, in which the boundary gradient condition is 
replaced by a (weaker) condition on the distance between two different level 
sets of the solution. More precisely, given the p-capacitary potential u in 
the annulus f2 \ K we ask that 

(1.1) dist(x, {u = i}) = a, vig an, 

where {u = 1} is the /-level set of u and I, A are two positive given constants. 

Condition (jl.ip was considered in |50j [5T| I52j . while studying diffu- 
sion problems with stationary isothermic surfaces; and independently by 
the third author in [60], where he considered a discrete version of Serrin's 
problem. This second approach was inspired by a parabolic free boundary 
model in finance in which the distance between the free boundary (that is 
the zero level set for the solution) and the location of a source is prescribed 
(c.f. [H [33 EH [20] for an overview). 

The novelty in the approach introduced here is that our problem with 
condition (jl.ip approximates the solution of the classical Bernoulli problem 
in the following sense: as A — > 0, any p-capacitary function u\ vanishing on 
dQ and satisfying (jl.ip with I = coX converges to a p-capacitary function 
u vanishing on <9f2 such that |Vu| = uj on dVt. In some sense, this new 
formulation can be understood as a discrete Bernoulli problem. 

Two different problems can be set up according to where the free bound- 
ary lies in the annulus: the exterior problem (Pe), considered in Section [3[ 
and the interior one (Pi), considered in Section [H It is clear that the geom- 
etry of the domains will play a crucial role in both problems. 

In this paper we deal with two different geometries: the convex and 
the star-shaped case. We address the questions of existence, multiplicity, 
uniqueness and regularity for the discrete Bernoulli problem, together with 
convergence to the usual Bernoulli problem in the convex setting. 

The main ideas in our paper follow the approach by sub- and super- 
solutions introduced by Beurling in [13]. For the classical Bernoulli problem 
in the convex case with p-Laplacian operator, questions of existence and 
uniqueness were addressed in [36j ETJ [38] . It is noteworthy that the use of the 
properties of the distance function (jl.ip avoids many of the technicalities: for 
example, higher regularity for solutions is an immediate consequence from 
existence theory, which is not the case for the classical Bernoulli problem. 

The most standard approach for the Bernoulli problem is either varia- 
tional or by singular perturbation. In the case p = 2, the seminal paper by 
Alt-Caffarelli [9] shows existence of a solution using variational techniques 
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without any geometric restriction on the set d{u > 0}. Moreover it is shown 
that the free boundary is flat, and consequently of class C 1,a up to a hyper- 
surface of Tl n ~ l -measure zero. For N = 2 they prove that the free boundary 
is globally analytic (see [32] ) . Recently in [19] it was shown that there exists 
no variational solution with singular free boundary for N = 3. The question 
remains open for 4 < N < 6, since in [27] existence of minimal cones is 
shown for N = 7. 

For p y£ 2 see [23] and [25]; in particular regularity is considered in 
[251 ES], where the authors extend the results of [9] (C ,a regularity near flat 
points). They prove that for N = 2 the free boundary is globally analytic 
for 2 — 5 < p < +oo with 5 an absolute constant. For N ^ 2 not many 
results are available: Petrosyan in |56] shows that for any N < k* — 1 and 
any 2 — 5 < p < 2 + 5 the variational solution has globally analytic free 
boundary, where k* is the critical dimension above which any variational 
solution has free boundary with singular points. 

If one adds a priori geometrical assumptions, the problem becomes more 
tractable. If the initial given domain is convex, for any 1 < p < +oo there 
exists a unique classical solution, and the free boundary is also convex and 
C 2 ' a for both interior and exterior case. This is done in \36\ [371 EH] using 
Beurling's method on sub- and super-solutions. 

Although the literature is not yet exhaustive, one would expect similar 
results if we relax the geometry of the domain to be star-shaped only (see, for 
instance [7J). The p-capacitary potential in a star-shaped annulus has star- 
shaped level sets, according to the rearrangement results of |41] ; see Section 
5 for a more detailed literature on star-shaped Bernoulli problem. However, 
two questions remain open: (i) Existence of a solution using Beurling's 
method, which, by uniqueness, would coincide with the variational one. (ii) 
The best regularity one expects for the free boundary. 

The paper is structured as follows: after recalling in Section [2] some 
preliminary results that will be used throughout the manuscript, we will 
consider the exterior and interior problem in the convex setting, see Sec- 
tions [3] and 131 respectively. Our problem will be then considered for star- 
shaped domains in Section [5j In Section [6] we further extend the problem 
to more general cases, such as non-constant distance constraint (including a 
constraint involving mean curvature), and two- and multi-phase problems. 
We also give the proof of a Brunn-Minkowski inequality. In the last section 
we show that the proposed model might give insights for new numerical al- 
gorithms that could more efficiently approximate solutions to the classical 
Bernoulli problem. Several open and tantalizing problems are also outlined. 



3 



2 Preliminaries 



For the readers' convenience we briefly recall in this section the main exis- 
tence and regularity properties for p-capacitary functions that will be used 
later. 

For 1 < p < oo, we define the p-Laplacian operator A p as 
A p u := div(|V-u| p-2 Vu), 1 < p < oo, 
and the equation A p u = in is understood in the weak sense as 

/ |Vn| p ~ 2 VuVv? = 0, V tp G WQ' p (n). 
Jn 

We first review some classical regularity properties: 

Theorem 2.1. Let U be a bounded domain inM. N andu be a (weak) solution 
of A p u = in U. Then: 

i. uecj£(U), (MW), 

ii. u satisfies the weak maximum principle and Hopf's boundary lemma 

m 

Additional Holder regularity up to the boundary can be shown through 
construction of conical barriers. With this objective in mind, we introduce 
some notations: we say that a domain U satisfies a uniform exterior cone 
condition if there exists a ro > such that for all x° £ dU the finite right 
circular cone V x o with vertex x° and opening ro satisfies UnV x o = {x } and 
V x o C U c . One may analogously define the uniform interior cone condition. 

Lemma 2.2 ([61]). Let U be a domain with exterior cone condition. Then 
the solution of the Dirichlet problem with Holder boundary data for the p- 
Laplacian in U is a Holder- continuous function up to dU with Holder norm 
uniformly bounded depending on ro, sup;y \u\ and Rq (radius of largest ball 
that fits inside U ). 

With some additional information on the regularity of dU one can obtain 
higher regularity for u up to dU, in particular, uniform gradient bounds. 

Lemma 2.3 (08]). Let U be a C 1 ' domain in R N . Then any solution to 
A p u = in U is C l ^(U) where j3 depends on p, N and a. 
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Finally, we say that a set U satisfies the uniform interior ball condition 
if there exists ro > such that for each x £ dU, there exists z x € U with 
the property that x € dB ro (z x ) n c?£7 and S ro (z x ) C U. 

In this paper we will be interested in the case that u is a p-capacitary 
potential in a convex annulus. Let K,Q be two domains in M. N , N > 2, 
satisfying if C O, f! bounded. In the sequel we say that un is the p- 
capacitary potential of the ring shaped region Q \ K if: 

( A p u n = in_ ft \ Z, 
(Pn) < «n = l onK, 

^ uq = on 

It is clear that a domain £l\K, with and K convex, automatically 
satisfies the uniform exterior cone condition, for some ro depending only on 
K. Thus, the Holder estimate from Lemma 12.21 will only depend on the 
(interior) cone condition for the domain K but not on Vt. We formulate this 
as a Corollary: 

Corollary 2.4. Let Vt, K be two bounded convex domains such that K C 
K C O, and denote by uq the p-capacitary potential of Q\K as given in 
equation (Pn)- Suppose that K satisfies the uniform interior cone condition 
for opening ro, and set do := dist(d£l,K). Then there is a constant M = 
M(ro,do) such that 

l|un|| c «»(n\F) < M - 

If, in addition, K is C 1 ' , one can obtain uniform gradient estimates up 
to the boundary: 

Corollary 2.5. Let Q,K be two bounded convex domains such that K dVL, 
and denote by uq the p-capacitary potential of Q\K as given in (Pn). Sup- 
pose that K satisfies the uniform interior ball condition for radius r$. Denote 
also do = mm x& dQ dist(x, K). Then there is a constant M = M(ro, do) such 
that 

|V«n|<M in n\K. 

Now we deal with the lower bounds for the gradient: 

Lemma 2.6 (Lemma 2 in 06]). Let £l,K be two bounded convex domains 
such that K C and denote by un the p-capacitary potential of £1\K as 
given in (Pn)- Suppose that £l\K satisfies the uniform interior ball condition 
for radius ro- There is a constant Mo such that 

[Vitfi| > M in n\K. 
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If K is uniformly convex, it can be easily shown by constructing a suitable 
barrier that the gradient of |V«q| near dK must be bounded from below 
by a positive constant. If we do not have this extra information on K, but 
we know the behavior near 0, we can translate this information back to 
K. First we need to introduce some notation. For two nested convex sets 
K C O, and for x € dK, we denote by T xa the supporting hyperplane at 
x with the normal a pointing away from K. Depending on the geometry of 
dK, T xa is not necessarily unique. Now for each x € dK there corresponds 
a point y x , not necessarily unique, on Oil. {z : a ■ (z — x) > 0} such that 
a ■ (Ux — x) = max a • (z — x), where the maximum has been taken over all 
z € dQ fl {z : a ■ (z — x) > 0}. Then, a crucial result is: 

Lemma 2.7 (Lemma 2.2. in |37J). Let K and O be two nested open convex 
domains, K C 0, and denote by uq the p-capacitary potential of £l\K , as 
given in problem (Pn)- Then 

limsup |Vii(z)| > limsup |Vu(z)| for all x € dK, 

zen\K zen\K 

where y x is the point indicated in the previous discussion. 

Next, a fundamental result on properties of level sets for p-capacitary 
functions uq in a convex ring: 

Proposition 2.8 ([46]). Let I € [0,1). Given a convex ring £l\K , its p- 
capacitary function uq described as above is real analytic and the level sets 
{uci > 1} are convex (with strictly positive principal curvatures) and analytic. 

Finally, a crucial technical result that allows to define the gradient of 
the solution at the boundary of a convex ring, and thus, give a notion of 
classical solutions for the usual Bernoulli problem: 

Proposition 2.9 (Theorem 1.3 in [35]). Let K and O be two nested open 
convex domains, K C £1, and denote by uq the p-capacitary potential of 
£l\K, as given in problem (Pn). Suppose also that dL( and 8Q, are C . 
Then us semi- continuous in Q\K, and non-tang entially continuous up 
to on U dK. 

Finally, we recall the definition of extremal points: 

Definition 2.10. Let £1 be a convex set: a point x € d£l is an extremal point 
if x cannot be written as a linear combination of the form x = tx\ + (1 — t)x% 
for t G (0, 1), xi,x-2 £ dQ. 



In the rest of the paper we denote by Yi the I level set of u: 
Yi := {x G Cl | u(x) = I}. 
We also write K d Cl when K C 0. 

3 The exterior problem in convex setting 

The first aim of this paper is to solve the following problem: given a convex 
open bounded set K C R , N >2, and / G (0, 1), A > positive constants, 
find a function ti and a convex open bounded domain Cl C M. N , Cl D K, 
solving 

{Apti = in Cl \ K, u = 1 in K, u = on <90, 
dist(x,r/) = A for all xe da 

We first show: 

Theorem 3.1. Problem (Pe) has a unique regular solution (u, CI) with Cl 
convex, dCl € C°° , and u G C 1,a (Cl\K) uniformly up to the boundary of Cl. 

The proof is inspired by previous works of one of the authors |36|, [37] 
via Beurling's method (see also the classical reference |13j). The method is 
based on constructions of super- and sub-solutions for (-Pe); hence, existence 
of a solution for (Pe) follows by the proof of the existence of a minimal 
super-solution which turns out to be a subsolution (or solution) as well. 



3.1 Proof of existence 

We first deal with the existence statement in Theorem 13.11 For each convex 
set Cl, let uq be the p-capacitary potential of Cl\K. Let us consider the class 

C := {Cl convex bounded open set in R N , Cl D K}, 

and define the following three subsets: 

A = {Cl G C | sup dist(x,T/) < A}, 

xedn 

A = {Cl g C | inf dist(x,r/) < A}, 
B = {Cl G C | inf dist(x,r;) > A}. 

Note that A is the set of subsolutions, *4o the set of strict subsolutions, 
while B is the set of supersolutions for problem (Pe)- Let us first check that 
they are non-empty sets : 
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Lemma 3.2. There exist a super solution and a strict subsolution for prob- 
lem (Pe)- 

Proof. By translation suppose that the origin of the cartesian system lies 
inside K. In order to construct a supersolution take < r < R big enough 
such that K C B r (0). The solution of the capacitor problem for the annulus 
Bji(0)\B r (0) is explicitly given by function 



log (f ) 

p—N p — N 

X p ~ 1 — ^ P — 1 

7" P- 1 — J? P-l 

It is clear that we can have dist({u = l},{u = 0}) as big as we wish by 
fixing r and taking R big enough. Consider also ur to be the p-capacitary 
potential of the set Br(0)\K. By the comparison principle, ur < u for all 
x G Bft(0)\B r (0). This immediately implies that 

(3.1) dist(x, {u R = /}) > A for all x € dB R (0), 

which shows that ur is a supersolution for our problem, i.e, Br(0) G 

Now we seek a strict subsolution. Let Br(0) and as above. Since 
K C Br(0), and is the p-capacitary potential in this annulus, we can 
apply Lemma 12.61 to bound the gradient of ur from below. More precisely, 
there exists a small neighborhood U of dK and some a > such that 

(3.2) |Vu fl (x)| > a > for all x€U\K. 

Choose the level set {ur = 1 — e} and set U e := {x : ur(x) > 1 — e} n U. 
Define now 

u R (x) - (1 -e) 
u e (x) := . 

The function u £ {x) is a p-harmonic function in the domain U e , u e {x) = in 
{ur{x) = 1 — e} and u £ (x) = 1 in K, so it is the p-capacitary potential of 
the set U e . We need to prove that u e G Ao, i.e., 

sup dist(x, {u £ = I}) < A. 

a:e{t%=0} 



This is easy to see since from the lower bound (|3.2p it follows 

dist(x, {u e = /}) < — < A, x G {u E = 0}, 
a 

by choosing e small enough. □ 
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Lemma 3.3. Let (us,£ls) be any super solution and (u s ,£l s ) any subsolution 
for problem (Pe)- Then £l s C £ls and u s < us- 

Proof. We use Lavrent'ev rescaling method. Assume by contradiction that 
there is a subsolution (u s ,Q s ) which is not smaller than a supersolution 
(wg, fig). Without loss of generality, assume that 6 K. Rescale the 
function u s by a parameter e < 1 small enough, i.e. u e s (x) = w s (|), so that 
u e s is a p-capacitary function in Q £ S \K E for £l e s C Og, K e C K. Moreover, for 
all x£dfl e s , 

dist(x, K = 0) < A £ < A. 

Consider the biggest £0 such that f^ C and dQl° n <9f2s 7^ 0. Note 
that we can take eo < 1 because £l s ^ tts- Then, by strong comparison 
principle u e s < us for each x G Q e s °\K. Let x° E dfll n dfis: it holds that 

A > A e > dist(x°, {u £ s = I}) > dist(x°, {u s = I}) > A, 

which is a contradiction. □ 

The next step in the application of Beurling's method is to show that 
the set B is closed under intersection. This is the preparation step in order 
to show that the limit of a decreasing sequence of elements in B still belongs 
to B. 

Lemma 3.4. Let ill and O2 be two elements of B. It holds that 

fiino 2 e B. 

Proof. Let itrhnf^ be the solution of the p-capacitary problem in the inter- 
section fi] n O2. By application of the maximum principle in each of the 
domains Vt\ and we can conclude that 

un 1 nn 2 (x) < mm{un 1 (x),un 1 (x)}, for each iGfiifl Cl 2 - 

Moreover, since Q\ n SI2 must contain K, the level set {un in n 2 = is weu 
defined. Let x° £ {uQ 1 nn 2 = ^} the point where the minimum distance 
between the sets d(Q± n ^2) and {un in n 2 = i s achieved. Then 

i = ^nf^ ) < min{«ni(aJ°),«n 2 (a; )}> 

which implies that the point x° is further away from both boundaries <90i 
and dfl2, in other words a distance greater than or equal to A. □ 
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Next lemma shows convergence of the sets O re . 

Lemma 3.5. Let fii D O2 D • • ■ 6e a decreasing sequence of convex domains 
in B, and suppose that the set Q, defined the as interior of the closure nflk, 
is convex, bounded, open and £1 D K. Then £1 € B. 

Proof. Let Uk be the p-capacitary potential of Qk\K. By construction, < 
< 1, and it is a decreasing sequence. Then we can easily show that on 
every compact subset of O := Interior (nS7fc), Uj~ converges in C 1,a norm to a 
p-harmonic function u. However, the lack of control of the gradients at the 
boundary does not allow to show higher regularity up to the boundary, but 
only C a . 

We first prove that u is precisely the solution to the p-capacitary problem 
in Q\K, denoted by uq. Clearly, u = 1 in K. We also need to check that 
u = outside f2. For that, let a; O. Since K is convex, £lk\K automatically 
satisfies the uniform exterior cone condition, and one may use Corollary 12.41 
to estimate 

(3.3) < u k (x) < dist(x,R N \ Sl fc ) Q sup [|t*fe llc«(0\J^) 

< Mdist{x,R N \n k ) a ->0, 

as k — > 00. This implies that u(x) = 0, as desired. 

Since the sequence of functions is decreasing, the set {u k = 1} converges 
to the set where {u = I}. However the boundary of the set may not 
converge to the set {u = 0}: that may happen if the gradient of u vanishes 
at some point of dtt so that the limit function flattens out, and in this case 
the distance property would be violated. As a consequence, we need to check 
that u(x) > for all x € f2. By the minimum principle, u cannot vanish in 
the interior, but it could happen that it vanishes at the boundary of some 
set C fi. Hence assume by contradiction that £1 ^ £1; consider a point 
x° € dCl, x° and let P € {u = 1} be the point where distance between x° 
and {u = 1} achieves its minimum. It can happen that x° = P but this is not 
a problem in the argument. Note that the ball B £ /2{ x °) is completely inside 
fife for all k (otherwise, the distance property for u k would be violated). 
Moreover u k {x) — > as k — >• 00 for all x € P> £ / 2 (x°) \ (B £ / 2 ( x °) H Cl). Cover 
the segment x°P by a finite number of balls of radius e/2, and define S = 
U^ i? e/ / 2 (x*). In the set S the Harnack inequality applied to the function 
u k implies 

/ = u(P) < Uk(P) < sup Uk(x) < C inf Uk(x), 
s s 
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where C depends only on S. But this is a contradiction since, as k — > oo, it 
holds 

I < C inf u k (x) -> 0. 

We have shown that O = 0, so that u > on fi. Then it = is the 
p-capacitary potential for the set Q\K. Moreover the set f2 belongs to B 
since 

(3.4) A < dist(<9fi n , {u n = I}) dist(90, {u = /}) as n — > +oo. 

Note here that the passage to the limit is justified because dQ n —¥ d£l and 
{u n = 1} — > {u = 1} in the Hausdorff distance sense. □ 

The previous two lemmas allow us to find a minimal set O in B, i.e., the 
solution u will be the minimal supersolution, that stays above the strict sub- 
solution we have found in Lemma [3 .21 More precisely we have the following 
corollary: 

Corollary 3.6. Assume that there exist two domains Q° £ *4o and Q 1 £ B, 
with n° C fl 1 . Set 

SS := {ft € B : fl° C fi}. 

Then there exists a domain Q in the class SS which is minimal for the 
inclusion. 

Proof. Let O be the intersection of all domains in the class SS and set Q to 
be the interior of the closure of 0, which is still convex. To prove £ B, we 
select a sequence of domains {Uk}^ = i in SS such that HUk = ^. Consider 
the sequence of domains defined by fii = U\ and £lk+i = ^fc H U^+i for 
> 1. Each 0^ is convex and because of Lemma 13.41 belongs to B. Lemma 
13.51 completes the proof of the corollary. □ 

In the following, we always denote by O the (non-trivial) minimal element 
in the class SS defined in the previous proposition. To finalize the proof 
of the theorem we need to check that this minimal set satisfies the distance 
property so it is indeed the solution of (Pe)- 

Proof of Theorem 3.1: 

We need to show that £1, the minimal set of B, and its corresponding p- 
capacitor potential uq satisfy the distance property: 

dist(x, {uq = I}) = X for all x £ dil. 
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The proof for the standard Bernoulli problem uses crucially the convexity 
assumption. Our problem may be solved through a different and much 
simpler argument, which may be extended to more complex geometries, for 
instance, the star-shaped case. By contradiction, suppose that there is a 
point x° G d£l such that 

dist({ii^ = I}, x°) = A + e, e > 0. 

Set £ := {uq > Note that C\K is a convex ring. Define £1 as the set 

:= {x G : dist(x, C) < A}, 

with u the p-capacitary potential of Q\K. By the hypothesis, f2 ^ and 
hence the comparison principle gives that u < u. One may estimate, for 
x G dCl, that 

dist(x, {u = I}) > dist(x, {u = I}) > A. 

Since <9f2 lies at fixed distance A from C, and £ is a convex bounded set, we 
can conclude that also O is convex. Hence belongs to the set B; but this 
contradicts the fact that f2 is the minimal set in B. 

Lemma 3.7. Assume that K is star-shaped. If problem (Pe) admits a 
solution, then it is unique. 

Proof. It is a consequence of Lemma [3.3| since a solution is at the same time 
a sub- and a super-solution. □ 

Remark 3.8. The uniqueness result together with the constructive construc- 
tion of Q as intersection of convex sets give that the solution O of problem 
(Pe) must be convex. 

We look now at regularity. It is clear that u G C°°(£l\K) and even 
analytic in the interior of Q\K, but what about the regularity of the free 
boundary? First note that {u = 1} is level set of an analytic function, hence 
it is analytic. Moreover all the normal curvatures at the level set are positive 
(Lemma l2.8p . The level set {u = 0} is at a fixed distance A from {u = 1} and 
its curvatures are of the form , where Hi are the curvatures of {u = I}; 
hence we can conclude [31] that {u = 0} is also C°° and analytic. 

The proof of Theorem 13.11 is thus completed. 
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3.2 Convergence to the Bernoulli problem 

Finally, we show that the discrete Bernoulli problem converges, at the limit, 
to the usual Bernoulli problem: 

Theorem 3.9. Let l n and \ n be two decreasing sequences of positive num- 
bers such that l n — > 0, X n — > as n — >■ +oo and l n = \ n uj. Let (u n ,£l n ) 
be the solution of the exterior problem (Pe) with constants I = l n and 
A = A n . As n — > oo, {£l n } is a decreasing sequence of sets converging in 
Hausdordff distance to some convex set Q and the function u n converges in 
C^(Q\K) nC a (Q\K) to the (unique) solution u of the classical Bernoulli- 
type problem 



(Pb) 



A p u = in f2 \ K, u = 1 in K, u = on d£l, 
|Vu| = cj, for all x G dQ. 



The boundary condition |Vu| = co is to be understood in the following sense: 
(3.5) liminf |Vu(y)| = limsup | Vii(ii) | = lj, for every x E d£l. 



Remark 3.10. It is shown in that the solution for the Bernoulli problem 
(Pb) exists, it is unique if K is bounded, and has C 2,a boundary dfl. 

Proof. Define the rescaled distance function 

d n {x) = - dist(x, {u n = Z„})/A„]. 

Clearly, when x E dVL n U {u Tl = l n } we have that d n (x) = u n (x). Define also 
S n '■= {0 < u n < l n }, a narrow ring. Since the level sets of u n are convex 
surfaces, the function d n is p-superharmonic in the set S n , i.e, 

-A p d n > 0. 

Since d n and u n have the same boundary values on dS n , comparison principle 
yields that u n < d n everywhere in S n . In particular, one obtains that for 
x E d£l n , and l m < l n , 

dist(x, {u n = l m }) > dist(x, {d n = l m }). 

On the other hand for any x E {0 < d n < l n } it holds \S/d n (x)\ = to. Since 

5; lm ^ri) 

dist(x, {d n = l m Y) > = — = X m for all m > n. 

sup|Va n | co 
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Thus, u n is a supersolution for every problem (Pe) with constants (A m , l m ), 
m > n, i.e., fi n G £> m where 

B m := {Q G C | inf dist(x, {uq = > A m }, 

for all m > n. Because of our construction of the solution for the exterior 
problem with constants (A m , l m ) as the minimal set in B m , this immediately 
gives that Q n contains all the Q m for m > n. Consequently we have shown 
that the sequence of {Cl n }^ =0 is decreasing and from maximum principle it 
follows that {u n } is a decreasing sequence. 

Since u n < d n on S n , using Hopf 's boundary lemma we arrive at |Vu n | < 
|V(f n | = uj on dQ n . This tells us that u n is a supersolution for the Bernoulli 
problem with constant uj in the sense indicated in |37j . 

As a consequence, the sequence of sets O n converges (for example in the 
Hausdorff measure) as n — > oo towards a convex domain Q D K. Moreover 
u n converges in C^°-norm to a p-harmonic function u: by Theorem 3.2 
in [37] we know that u solves the p-capacitary problem in \ K and is a 
supersolution, i.e., |Vu| < w for all x in dQ. 

The last step consists of showing that (u, Q) solves the Bernoulli problem 
(Pb) with constant oo in the sense indicated in (j3.5|) . Use again comparison 
to the distance function, we obtain that |Vu n | > uj on the level set {u n = l n }. 
But Lemma 12.71 gives that the inequality | Vu„ | > u is also true for every 
point in {l n < u < 1}. Since we have uniform convergence of the gradients 
away from the boundary and the sequence of sets O n — > Q, is decreasing, we 
are able to pass to the limit and thus |Vn(y)| > uj for every point y € Q\K. 
The proof of the Theorem is completed. □ 



4 The interior problem in convex setting 

We consider now the following problem: given a convex open bounded set 
Q C 1^ for N > 2 and I G (0, 1), A > 0, find a function u and a convex 
open bounded domain K C f2 such that 

{A„u = in J7 \ K, u = in K, u = 1 in dtt, 
dist(x, T[) = A, for all x G dK. 

Our first main result for the interior case is the following. 

Theorem 4.1. There exists a constant \n,max that depends only on I and 
Q such that for any A < \n,max problem (Pj) has a smooth solution (u,K) 
with K a convex set and dK G C ' 1 . 
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For each convex set K, one can define its associated p-capacitary poten- 
tial uk as the solution of the following problem: 



CP*) 



' A p u K = in O \ K, 
uk = in K, 
uk = 1 on d£l. 



Note that, without risk of confusion, we have switched the notation from 
problem (Pq) in the exterior case. 

The proof of Theorem 14.11 goes, as in the exterior case, by considering 
the maximal element in the class: 

B\ = {K convex bounded open set, 

( 41 ) n D K | inf dist(x, {u K = I}) > A}. 

x£dK 

In the following we show that, fixed I € (0, 1), the set B is non-empty if the 
value of A is smaller than a certain critical value, which will be denoted by 
^n,max- We also show that problem (Pi) does not necessarily have a unique 
solution for certain values of A. 



4.1 A minimal supersolution and the Bernoulli constant 

As a motivation, consider the simple case that f2 is a ball: such particular 
example shows, as in the case of the classical Bernoulli problem, the existence 
of a constant Xn,max such that no solution exists for A > \n,max, one or two 
solutions can occur for A < \n,max and only one solution exists if A = \n,max- 
Using Beurling's terminology, we call the unique solution that corresponds to 
the value A = \n,max of parabolic type, the smallest solution for A < Xn,max 
of hyperbolic type and the biggest solution for A < \n,max of elliptic type. 
The elliptic solutions form a decreasing family, while the hyperbolic solutions 
form a increasing family of solutions. 

Lemma 4.2. Let Q = B R (0) CR N be a ball with radius R. Fixed I € (0, 1) 
and p > 1, we have that: 

i- If p < N , there exists a constant X max depending on I and R such 
that problem (Pi) has a unique solution for A = \ m ax> two solutions if 
< A < Xmax and no solutions if A > X max . 

ii. If p ^ -/V, then there exist constants A ma ^ and X m i n depending on I, p 
and R such that problem (Pi) has a unique solution for X G (0, X m i n ) U 
{X max }, two solutions if X £ (Xmin, Xmax) o,nd no solutions if X > 

Xmax ■ 
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Proof. It has been shown in [57] that all the solutions of the p-capacitary 
problem with £1 = -Br(O) are balls centered in the origin. We deal first with 
the case N = p. Given any r with < r < R, the unique solution u r of the 
p-capacitary problem with f2 = -Br(O) and K = B r (0) reads as follows: 

log(.R) - log(r) 

It holds 

u r (|x|) = I if and only if \x\ = R l r l ~ l , 
and consequently 

A(r) := dist(x, {u r = l}) = r ((f)'- l) for all x 6 dB r (0). 

We note here that A(r) + r < for r S (0, i?) so u r is an admissible solution 
for the p-capacitary problem in the annulus. 

With elementary computations one can check that A(r) reaches its max- 
imum value at 

i (l-l) 1 / 1 
(4.2) r max :=R(l-l)i with A(r max .) = _ ^ =: A maa; . 

We may easily compute A(0) = 0, A(R) = 0, and A(r) is strictly increasing 
for < r < r max , strictly decreasing for r max < r < R. This implies that for 
any A G (0, \ m ax) there exist two values rx,i < r max < r\ t 2 < R such that 
both solutions Ui, i = 1,2 of the p-capacitary problems in -Br(O) \ B rxi (0) 
satisfy dist(x,{uj = I}) = A for all x € dB Ti (0) for i = 1,2. Moreover, 
there exists a unique solution for X ma x, and no solutions in the remaining 
case. This completes the proof of the lemma for the case p = N. 

On the other hand, when p ^ N the p-capacitary potential is given by 

p-JV p—N 

r p- 1 — \x\ p- 1 

U r (\X\) = ^Tjy p-JV ; 

f p— i — R v— 1 
and, keeping the same notation as above, 



p-JV p— JV 

A(r) = (l-l)rp- 1 ilflp- 1 



p-i 

p-JV 

— r. 



But again, A(r) + r < R, so u r is indeed an admissible solution in the 
annulus. 
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Next, A(r) reaches its maximum value at 



R 



N-p 



p-N 



Note that A is ctn increasing function in the interval (0, T max ), decreasing 

at (r max ,R), A(R) = and A(0) = R > if p > N and A(0) = for 
p < N. 

We set Xmin '■= A(0) and 
(4.3) 



/ \ p-N 



i ) p - N _ L 

P-l I N-p 



".1-(1-Z)*T=T/ ^(l_J|)7r=T 

This completes the proof of the lemma. □ 

Now we go back to the problem (Pi) for a general convex bounded set f2 
in M. N and a fixed / G (0, 1) (although we may not write it, the dependence 
in I will be always implicit). Consider the set of supersolutions as defined 
in (I4.ip . Since B\ x C B\ 2 for Ai > A2, we can define the Bernoulli constant 
for our problem as 

(4.4) \n,max ■= sup{A > I B\ is not empty}. 

A possible bound from below for \£i,max is given by the next lemma. 
Lemma 4.3. Let I G (0, 1) and £1 be any open bounded convex set, then 

Af2, max — ^maxt 



where X ma x is defined in (I4.2p if p = N and (|4.3p if p ^ N, and R is the 
radius of the largest ball inscribed in £1 . 

Proof. According to Theorem 14.11 and the definition in (|4.4j) . problem (Pi) is 
solvable for given constants A > 0, / G (0, 1), provided that B\ 7^ 0. Follow- 
ing the discussion above, for any A < A(r max ) there exists a unique constant 
f\,2-, with r max < ta,2 < R such that the solution of the p-capacitary problem 
u rx 2 in -Br(O) \ B rx 2 (0) satisfies the distance property: 

dist(x, {u rx 2 = /}) = A for all x G dB rx 2 (0). 

Note that if A = D(r max ) then r max = r A>2 . 
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For any A < D(r max ) consider now the solution u of the p-capacitary 
problem in Q \ B rx 2 (0). By comparison principle applied to the functions 
u rx 2 and u, the (minimum) distance between the sets {u = 1} and dB rx 2 (0) 
is greater or equal A. This shows that u is a supersolution for (Pi), as 



Corollary 4.4. Fix I € (0, 1) and Q a convex bounded domain in M. . Then, 
for each A < \n,max, there exists a supersolution (u, K) such that K belongs 
to the class B\. 

Proof. Indeed, we have just seen that if xq £ Q and R the radius of the 
largest ball centered at xq inscribed in fi, then for A < \n,max the solution 
u constructed in Lemma 14.21 is a supersolution of the problem and indeed 
its zero set K := B rx 2 (0) belongs to the class B\. □ 

Remark 4.5. Because of the distance property, for any K 6 B\, we have 
that 

dist(x, dQ) > A for any x € dK, 

which shows that any element in the class B\ is strictly contained in f2 and 
thus, non- degenerate. 

We have just shown that for any A < Xn,max the set B\ is not empty. 
Obviously it holds that B\ 2 C B\ x if Ai < A2. We first define the maximal 
set as 



and will prove that KL\ satisfies the distance problem (Pj) with constants 
(A, Here ^(X) denotes the closed convex hull of the set X. 

It is an interesting open question to see if for general convex domains 
ft the set of all solutions of problem (Pj) follows the same structure as the 
ball case, i.e. parabolic, elliptic and hyperbolic case; in particular to prove 
that the family of maximal convex solutions defined above is an elliptic 
(i.e. decreasing and continuous) family of solutions and that there exists 
a unique solution corresponding to the value A = Xn )ma x- These questions 
have been already addressed for the classical Bernoulli problem and proven 
to be true for the linear case p = 2, [31 [21] and for p ^ 2 more recently in [15] . 

Before we prove Theorem I4.lt let us show a couple of preparatory lem- 
mas: 



stated. 



□ 



(4.5) 
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Lemma 4.6. Let K\ and K2 he two elements of B\. It holds that 



K* € B x , 

where K* is the convex hull of the set K\ U K2. 

Proof. Let u* be the solution of the p-capacitary problem in \ K* . Since 
the set K* is strictly contained in £1, the level set {u* = 1} is well defined. 
Let x° 6 {u* = 1} and x 1 6 dK* be the points where the minimum distance 
between the sets dK* and {u* = 1} is attained. We would like to show that 
disttx^x 1 ) > A. 

First, by application of the maximum principle we can quickly conclude 
that 

u*(x) < min{uK 1 (x),UK 2 ( x )}i f° r each x G Q, 
which, in particular, implies that 

I = u*(x°) < min-jwifj (x°), uk 2 i x °)}- 

Because, by initial hypothesis, uk x and uk 2 are supersolutions of the prob- 
lem (Pj) in the sense of (|4.ip . one knows right away that every point in 
{u* = 1} is far from both boundaries dK\ and 8X2 a distance greater or 
equal than A, and in particular the same is true for x°. We have just shown 
that if we set K := K\ U K2, then 

(4.6) dist(x°,dK) > A. 

On the other hand, we have constructed K* as the convex hull of K. 
If x 1 is an extremal point of K*, then it belongs also to dK which implies 
that dist^ ,^ 1 ) > A by flM]). If x 1 is not an extremal point of K*, it can 
be written as a linear combination of points {y«}" =1 that are extremal for 
K. For such points we have already shown that dist(yj, {u* = I}) > A. But 
because the set {«* = /} is the boundary of a convex set by Proposition 12.81 
one automatically obtains that dist(x 1 ,x°) > A too. □ 

Lemma 4.7. Let K\ C K2 C ... he a increasing sequence of convex domains 
in B\, and suppose that the set K, defined as interior of the closure L)K n , 
is convex, bounded, open and f2 D K. Then K G B\. 

Proof. Let u n be the p-capacitary potential of Q\K n . By construction, 
< u n < 1, and {u n } is a decreasing sequence. Then we can easily show that 
on every compact subset of O, u n converges in C 1,a norm to a p-harmonic 
function u. We first prove that u is precisely the solution to the p-capacitary 
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problem in £l\K, denoted by uk- Clearly, u = 1 on dtt. We also need to 
check that u = in K. For that, let x G K \ K n . Since K is convex, £l\K n 
automatically satisfies the uniform exterior cone condition, and one may use 
Corollary 12.41 to get a uniform (in n) C a estimate up to the boundary, hence 

< u n (x) < dist(x, K n ) a sup\\u n \\ Ca (ji\ K ^ < Mdist{x,K n ) a -> 

as n — > co. This implies that u(x) = for all x G K, as desired. 

Note that, because we have a decreasing sequence of functions, the set 
{u n = 1} converges to the set where {u = I}. However, the boundary of the 
set K n may not converge to the set {u = 0}, and this may happen if the 
gradient of u n grows in some neighborhood of dK. Hence we need to check 
that u(x) > for all x € £1 \ K. Assume that K C K and u = in K . This 
is done again by a Harnack chain argument: let x° G dK and P € {u = 1} 
be two points such that 

dist(x°, P) < A — e for some e > 0. 

Since K n satisfies the exterior ball condition, the ball B £ / 2 (x°) is completely 
outside K n for all n. Moreover, u n (x n ) — > for all x G B £ /2(x°) as n /oo. 

Cover the segment x°P by a finite number of balls of radius e/2, and 
define S = ^f = iP> e / 2 (xi). In the set 5, Harnack inequality applied to the 
function u n implies 

/ = u(P) < u n (P) < supu n (x) < C m£u n (x), 
s s 

where C depends only on S. But this is a contradiction, since as n / oo, 

I < CmfuJx) -> 0. 
S 

We have shown that, necessarily, K = K, so that u > on \ K. Then 
u = uk is the p-capacitor potential for the set Q\K. Moreover K belongs 
to B\; the proof follows the same argument as in (|3.4[) with Cl n replaced by 
K n . □ 

Now we check that the maximal set found in the previous lemma satisfies 
the distance property. Such maximal set provides the solution for problem 
(Pi), and thus the proof of Theorem 14.11 is almost complete. 

Lemma 4.8. Assume that B\ is non-empty. Let K be the maximal set of 
B\ found in the previous lemma and let uk be the corresponding solution of 
the p-capacitor problem for Q\K. It holds 

dist(x, {uk = I}) = A for all x G dK. 
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Proof. We know that both the sets dK and {u = 1} are convex (the sec- 
ond one is convex by Lemma l2.8p . Then one just needs to understand the 
distance between two convex sets. By construction, we know that for all 
x e dK. 

(4.7) d\st({u K = I}, x) > A. 

Suppose that there is a point x° € dK such that 

dist({u K = I}, x°) = \ + e. 

Then one can consider the convex set K = {x £ £1 : dist(x, {uk = I})} = A. 

Let u be the p-capacitary solution in K; clearly u < uk by the strong 
maximum principle. Then, if x € dK, 

dist(x, {u = /}) > dist(x, {uk = I}) > A, 

which shows that K is in the class B\, which contradicts the maximality of 
K. □ 

Conclusion of the proof of Theorem \4-l\ We have already seen that for 
any A < \n,max, the set B\ is nonempty and thus there exists a solution for 
(Pi). It remains to check that this is so for A = \n,max- 

Let Afc be a sequence converging to \n,max from below. We have just 
shown that for each A^, there exists a convex solution K^. One can subtract 
a subsequence, still denoted by Kk, which converges to a convex domain K in 
the Hausdorff topology. Then, by Arzela-Ascoli, the p-capacitary potential 
Uk of £l\Kk converges to u, the p-capacitary potential of £l\K. We have also 
shown that the distance from the l-level set of Uk is at a distance Xf~, then 
passing to the limit we also see that the same holds for u (similar argument 
as before). 

Remark 4.9. In contrast to the exterior problem studied before, in the in- 
terior discrete Bernoulli problem the free boundary may have corners; con- 
sequently dK is only C 0,1 and the solution is C°' a up to the boundary. 

4.2 Convergence to the Bernoulli problem 

Before we state the main theorem of the section, let us recall what is known 
for the interior Bernoulli problem: 
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Theorem 4.10 (|_36j). Given a convex domain £1, there exists a constant ujq 
depending on Q, and p such that for all uj > ujq the Bernoulli-type problem 



has a solution (v,D) (not necessarily unique). Moreover D is convex. The 
boundary condition |V^| = uj is again understood in the sense of (j3.5[) . 

Remark 4.11. We remind the reader that the solution (v,D) given in the 
previous theorem is the one constructed as the extremal supersolution for 
problem (Pbj), although there may be others. The set of super solutions for 
each fixed uj is defined in [36] as the class of Lipschitz functions in Q such 
that v = 1 on dQ, A p v < in {v > 0} D and \Vv\ < uj on d{v > 0} n 

Now, the main theorem of this section is: 

Theorem 4.12. Let Q be a bounded convex domain in M. N . Let uj > ojq, 
where ujq is the Bernoulli constant for the domain O given in Theorem \4-10 
Fix two decreasing sequences {l n } and {A n } of positive numbers such that 
l n — > 0, A n — > as n — >• +oo and l n = \ n ui. Let also ^Q max be the 
corresponding quantity (|4.4p to each l n . Then: 

• There exists a solution (v,D) for the interior Bernoulli problem (Pbi) 
with constant uj. 

• There exists a sequence of solutions (u n ,K n ) of problem (Pj) for each 



• The sequence of convex sets {K n } is decreasing and as n — > oo con- 
verges in Hausdorff distance to the convex set D. 

• In addition, the function u n converges in C,'^(fl\D) r\C a (Q\D) to v. 
The proof of the theorem will be a consequence of the following claims: 

Claim 1: The existence of a solution for the Bernoulli problem with 
constant u > ujq is the statement of Theorem 14.101 

Claim 2: For each n, there exists a solution (u n ,K n ) to the interior 
problem (Pj) with constants A = A n and I = l n . 

The proof of this claim uses the following lemma on the relation between 
the solvability of the distance problem and the Bernoulli one: 





A p v = in Q \ D, v = in D, v = 1 in d£l, 
\Vv\ = uj for all x 6 dD, 
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Lemma 4.13. Given a convex bounded domain £1 and I € (0, 1), the Bernoulli 
constant defined in (|4.4p for the distance problem satisfies 



(4.8) ^tt,max > — j 

where ujq is the Bernoulli constant from Theorem \4-10[ 

Proof. As in the proof of Lemma [4. 3\ it is enough to show that any solution 
(v,D) for the Bernoulli problem (Pb) with constant uq is also a supersolu- 
tion for (Pj) for A < Indeed, the Bernoulli condition for v tells us that 

|Vw| = <^o on 9D, which implies that \Vv\ < ojq in Q\D because of Lemma 
O Then, if x G dD, 

(4.9) dist(x, {u = /}) > l — — = — > A, 

K ' sup|Vv| w - 

which implies D £ B\. □ 
Using the lemma for the domain and the level l n we may estimate 

\ n > ^ n > ^ n — X 

A !],max — — — ^n, 
UlQ 0J 

which tells us that the distance problem (Pj) with constants (X n ,l n ) may 
be solved (Theorem I4.ip . This proves Claim 2. 



Claim 3: K n is a decreasing sequence of sets. In particular, we have that 
u = sup n u n and K = n n K n , where (u n ,K n ) is the solution to the distance 
problem (Pj) with constants A = A n , I = l n . 

Fix m < n, so that l n < Z m and A n < A m . Let us check first that u n 
belongs to the set of supersolutions B m , defined as 

B m := {K C | inf dist(x, {u K = l m }) > A m }. 

For that, consider the distance function d n between two level sets of u n , 
rescaled by a factor lj, so that d n and u n have the same boundary values in 
the set S n := {0 < u n < l n }. More precisely, 

(4.10) d n (x) = l n [l- dist(x, {u n = l n })/\ n ]. 

Since the level set {u n = l n } is smooth, the distance function is smooth 
(a reference for the regularity of distance functions may be found in j31j). 
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Note that in order to have d n well defined outside S n , we need to take the 
distance function to the level set dist(x, {/„ = u n }) with a change of sign 
once we pass the level set. 

Next, since the level sets of u n are convex surfaces (see Lemma l2.8p . the 
distance function defined in the set S n is p-subharmonic, i.e., 

A p d n > 0. 

Moreover, since d n and u n share the same boundary values on both bound- 
aries of the set S n , we can use Hopf's lemma on the smooth curve {u n = l n }, 
thus obtaining that 

uj = |Vd n | > |Vu n |, V x € {u n = l n }. 

Define now the set A n := {l n < u n < 1}. From Hopf's lemma we can 
deduce that the function f n '-=d n — u n is strictly positive in a annulus M n 
(a neighborhood of {u n = l n }) inside A n . Suppose, by contradiction, that 
M n is strictly contained in the neighborhood A n \ consider f n in the domain 
M n := S n U M n : it holds that f n = on dK n , f n > in {0 < u n < /.„}, 
f n = on {u n = l n }, f n > in M n and f n = on dM n . But this is 
a contradiction since f n is super-harmonic and should reach its maximum 
value at the boundary of M n . Hence M n = A n - Then in the set A n one has 
that u n < d n - As a consequence, for m < n, i.e., l m > l n , and x 6 dK n = 
{u n = 0}, 

dist(x, {u n = lm}) > dist(x, {d n = l m }) = X m , 

and u n € B m as claimed. In particular, since u m is constructed as the 
extremal supersolution in the set B m , we have that u m < u n , which in par- 
ticular implies that K n C K m - This is, K n — > K is a decreasing sequence of 
sets, and the claim is proved. 

Claim 4-' Let K := n n K n . Then K D D and is convex. In particular, 
it is not degenerate. 

Following the steps of Lemma [4.13[ one can show that the solution (v, D) 
of problem (-Pb 7 ) is a supersolution to the (A n ,/ n ) distance problem. Hence 
K n H D for any n 6 N. 

Claim 5: There exists a function u such that u n — > u in the space 

cl£(n\K)nc a (n\K). 

Claim 6: In addition, u = on dK and dK G C 0,a , so that u is the 
p-capacitary potential in the annulus Q\K. 

Claims 5 and 6 follow similarly as the exterior case as a consequence of 
uniform C a estimates up to the boundary for a convex annulus £l\K, so we 
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will not make further comment. 

Now we remark that the the regularity of the boundary dK is slightly 
worse than in the exterior case (see Remark 14. 9p . indeed it is only C ' 1 , and 
we may not have gradient bounds up to the boundary. However, we still 
have from the proof of Claim 3 that in the set A n := {l n < u n < 1} it holds 
that 

(4.11) \Vu n \<uj. 

Now we may pass to the limit n — > oo. It is easy to see that |Vu n | converges 
uniformly to |Vu| on compact sets inside Q\K. As a consequence of (I4.11j) 
we get the following claim: 

Claim 7: |Vu| <won dK. From Lemma [2.71 we also get that |Vu| < lo 
in tt\K. 

To complete the proof of Theorem 14.121 it suffices to show that D = K 
and u = v. First, in Claim 4 we have shown that D C K. On the other 
hand, in Claim 7 we have shown that |Vu| <won dK. This in particular 
implies that (u, K) is a supersolution for the Bernoulli problem (Pb 7 ) in the 
sense of Remark 14. Ill As a consequence, one gets that K C D (see |36j). as 
desired. 



Remark 4.14. If£l is a ball, relation (j4.8|) can be written explicitly. Indeed 
for \ m ax in Lemma \4-%\ via the limit 

i { e/R if N = p, 

lim = < i /„_i\ fzr^ 

one recovers the usual Bernoulli constant for the ball calculated in 136\/ . 



5 The star-shaped case 

A domain U C M. N is called star-shaped with respect to the point Xq 6 U if 
for all x € U the segment joining x and xq is contained in U. Moreover the 
set U is called star-shaped with respect to the ball B C U if U is star-shaped 
with respect to every point xq G B. 
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Let U be a star-shaped domain with respect to the point xq. For 9 G 
S^ -1 , where S^ -1 is the unit sphere, set 

$(0) = sup{p > : x + p9 g U}. 

Then we can write 

U = {x = p9-.ee § N ~\0 <p< (f)(9)}. 
It is well known j!7j that: 

Lemma 5.1. Let U be a bounded domain in M. N which is star-shaped with 
respect to the point xq G U. Then it is star-shaped with respect to a ball 
centered at xq if and only if the function <j) satisfies the Lipschitz condition 
on S N ~ 1 . 

We call such domains star-shaped Lipschitz. It is clear that star-shaped 
Lipschitz sets satisfy both the interior and exterior cone condition. One may 
also consider star-shaped Lipschitz rings: they are of the form U = Q\K, 
where f2 and K are star-shaped Lipschitz with respect to the same center. 

Given a bounded, star-shaped Lipschitz domain K C M N , we consider 
the exterior distance problem 




A p u = in £1 \ K, u = 1 in K, u = on dfl, 
dist (x,Ti) = A for all x € d£l, 



for some function u and a star-shaped domain £1 C M. , with the same center 
as K. Our main result in this section is stated in the following theorem: 

Theorem 5.2. Problem (P*) has a unique solution (u, Q) with dil G C°' a 
and C 1 ' 1 from inside. 

Before we give the proof of the theorem, let us briefly recall what is known 
for the (exterior) classical Bernoulli problem (Pb) when K star-shaped with 
respect to a ball. For the case p = 2 Kawohl proved in [41] existence and 
uniqueness of a (variational) solution via rearrangement method. In addition 
he proved that such solution has star-shaped Lipschitz level sets. Also Acker- 
Meyer [7] showed that, still for p = 2, the solution is regular and that 
d£l G C°° is also starlike with respect to the same ball. 

Now we are ready for the proof of Theorem 15.21 it follows the lines of 
the problem in the convex case, but some additional ingredients are needed 
to compensate for the lack of convexity. First, note that uniqueness was 
already shown in Lemma 13.71 
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Note that the regularity results of Theorem 12.11 and Lemma [22] still hold 
since they only depend on the Lipschitz constant of the boundary of the 
domain. Moreover, Lemma 12.61 is true for star-shaped rings as it was shown 
in [17]: indeed, a solution of the p-Laplacian equation is smooth, C 1,a , and 
satisfies a type of Harnack inequality, from which a bound from below for 
the gradient follows. Moreover, we have the following proposition: 

Proposition 5.3 ([41]). Let I G [0,1). Given a (Lipschitz) star-shaped 
ring tt\K and its p-capacitary function uq described as above, the level sets 
{uq > 1} are (Lipschitz) star-shaped. 

Since, by hypothesis, the initial domain K is star-shaped Lipschitz with 
respect to some point, then there exists some small ball B$ centered at that 
point such that K is star-shaped with respect to that ball. Consider the 
class 

C* := {O open, bounded, star-shaped with respect to B$, 0, D K}, 
and define the set of sub- and super-solutions as follows: 
A = {0, G C* | sup dist(x,r^) < A}, 

xedQ 

Ao = {n E | inf dist(x,rj) < A}, 
B = {ft G | inf dist(x,r z ) > A}. 

x£dQ 

The construction of super- and sub-solutions in Lemma 13.21 is valid also 
for K star-shaped Lipschitz domains; this implies that the sets Ao and B are 
nonempty. Lemma 13.41 holds as well. Consider now a set f^o £ Ao and set 
C M. N to be the intersection of all domains in the class B containing f^o- 
By construction, f2 is star-shaped with respect to the same ball B$. Let u to 
be the p-capacitary potential of Q\K. Then one may conclude as in Lemma 
13.51 that f2 G B; the only tricky point is to check that there exists a uniform 
Holder estimate in equation (|3,3p . But this is true because C Q -estimates up 
to the boundary from Corollary 12.41 only depend on the Lipschitz constant 
of the boundary. Moreover, we know by construction that 

inf dist(x,r ; ) > A. 

This implies that the equivalent of Lemma 13.51 and Corollary 13.61 for star- 
shaped domains are proven. 

Finally, one needs to check that the function u is indeed a solution for 
the distance problem (-P*), i.e. 

(5.1) dist(x, Ti) = A for all x G d£l. 
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This follows from a simple argument: let f2 be the extremal set as con- 
structed above, and set C := {u > I}. Note that C is a star-shaped 
annulus. By contradiction, assume that there exists xq G dQ such that 
dist(xo,5£) > A. Define Cl as 

Q := {x £ Q : dist(x, C) < A}, 

and let u be the p-capacitary potential of £1\K. By hypothesis, $7 ^ ft. 
Comparison principle implies that u <u and consequently 

dist(x, {u = I}) > dist(x, {u = I}) > A, for x G dCl. 

Since d£l has a fixed distance A from Ti, and T/ is a star-shaped surface, 
we conjecture that is also star-shaped with respect to the ball Bg. We 
proceed by contradiction: suppose that Cl is not star-shaped with respect 
to a point in B$. Without loss of generality we can assume this point to be 
the origin x = 0. Note that both sets C := {/ < u < 1} and £1 are smooth. 
Then there exists a point z G dQ such that the tangent plane II to dQ at z 
passes through the origin. 

By construction the point z lies at a distance A from Tf, in particular 
there exists xq G Ti where this distance is attained. Since all the sets here 
are smooth, we know that the line joining xq to z is normal to the plane 
II. Then one may just work with the intersection of Vt with the plane IIi, 
where IIi is the plane that contains the origin, the point z and the normal 
vector v (see Figure HJ. 

Consider now the line s passing through the points and z and the line 
t passing through and xq. Since the set C is star-shaped with respect to 0, 
the curve Ti n IIi has to cross the line t exactly at xo, and after crossing, it 
lies on the half-plane bounded by t that does not contain the point z. Hence 
there exist some points in O n IIi which are on the half-plane bounded by 
s (which does not contain xo) that lie at a distance strictly greater than A 
from T[. But this contradicts the definition of 0. 

Hence is star-shaped and belongs to B; by the minimality of Q if 
follows that O = f2 and (|5.ip is proven. 

For the regularity, just note that the distance property implies that <9f2 
is C ' from inside, but further regularity is unknown. The proof of Theorem 
15.21 is thus completed. 
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Figure 1: Star-shaped domain 



6 Generalizations 

6.1 Non-constant distance property - convex case 

The next question one may ask is what happens when we replace the con- 
stant A in the distance property by a non-constant function A(x). For the 
Bernoulli problem, this problem was considered in [35] • We concentrate on 
the exterior problem here; the interior problem follows the same line with 
appropriate modifications. 

Theorem 6.1. Let K C R N , N > 2, be a convex open bounded domain, 
I G (0, 1) and A(x) a continuous function satisfying < Co < A(x) < c\ in 
M. N . Suppose moreover that \{x) is locally concave in M. N \K (i.e., concave 
on each line segment contained in the setM. N \K). 

Then there exist a function u G and a convex open bounded 

domain C M. , D K, d£l € C 2 which solve the following problem: 



is strictly increasing for all x £ R for some zq G K, then £1 is unique. 

Proof. The proof of the above theorem is analogous to one of Theorem 13.11 
with small modifications that we indicate below. First, the sets A, B and 




A p u = in $7 \ K, u = 1 in K, u = on dQ, 
dist(x,Ti) = A(x) for allx G dQ. 



Moreover if the function 



(6.1) 



t G [0, 1] ^ t\((x - z )/t) 
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Ao are defined as in Section [3l Super and subsolutions are constructed as 
in Lemma 13, 21 by choosing, for the supersolution, R big enough such that 

dist(x, {ur = I}) > c\ 

and, for the subsolution, e small enough such that 

el 

— < co. 

a 

Next, Lemma 13.41 and 1 3.51 follow similarly. The existence of a neighborhood 
of radius e contained in each Qj~ for all k where one can apply Harnack 
inequality is a consequence of the fact that X(x) is a continuous function. 

Corollary 13.61 holds as well for A = X(x) and the proof follows exactly 
the one for A = const. Hence we obtain (O, u) minimal in the class B. is 
convex by construction, and thus, d£l is Lipschitz. 

It remains to check that this solution satisfies the distance property, i.e., 

dist(x, {u = I}) = X(x) for all x G dQ. 

Contrary to the constant case, convexity for is a crucial hypothesis. By 
contradiction, assume that the above is not true, i.e. there exists a point 
xq € 80, such that 

dist(xo, {u = I}) > X(xq). 

By continuity one can find a neighborhood V 3 xq such that dist(x, {u = 
I}) > X(x) for all x £ V. Let Eq be the set of extremal points of Q. If 
€ Eq, then one may find a convex set Q such that Q\V £ f2 $1 Q and 
xo £ Q. Let u be the p-capacitary function of £1. Comparison principle 
implies that u < u; then for all x € dVL, 

dist(x, {u = I}) > dist(x, {u = /}) > X(x). 

We have just shown that f2 G B. Contradiction with the minimality of Q. 

If xo Eq, then xq can be written as a linear combination of extremal 
points, i.e., 

n n 

(6.2) xo = ^^t n x n , where ^ tj = 1, Xj € E n ,i = 1, . . . , n. 
i=i i=l 

By hypothesis dist(xj, {u = I}) = A(x«) for every i = 1, . . . , n and 

dist(xo, {u = I}) < } j Udist(xi, {u = I}) = } UX(xi) < A(xo), 
i=i i=i 
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where we have used the local concavity assumption on A for the last inequal- 
ity. 

Next, to show that that the solution is unique we follow the same rescal- 
ing method as in Lemma 13.71 Suppose, by contradiction, that (P' E ) has 
two solutions (ui,f2i), (^2,^2)- Without loss of generality assume that in 
hypothesis (|6.1|) we have that zq = € K. Rescale u e 2 {x) = U2(x/e). Then 
Sl| C r^i and one may consider the biggest e(< 1) such that f2| touches f2i 
from inside. Let Xq € (9f2| H and xo = Xq/e. It holds 
(6.3) 

dist(xQ, {u E 2 = I}) = edist(xo, {^2 = I}) = eA(xo) = e\{x e Q /e) < A(xq), 

thanks to (16. ip for the last inequality. On the other hand, by comparison 
principle, we easily have that u\ < u\\ in particular, 

(6.4) dist(^, {u e 2 = I}) > dist(x§, {m = I}) = \{x e ). 

Then equations (I6.3P and (16.4p give the desired contradiction. 

We look now at the regularity: since is convex the surface {u = 1} is 
analytic by Lemma 12.81 Then dtt is the set of points at distance A(x) from 
this level set, that will have principal curvatures given by the formula 

1 + \{x)m' 

If A is continuous, then d£l will have continuous principal curvatures, which 
implies that <9$7 G C 2 . Higher regularity on A will imply higher regularity 
for dO,. The proof of the theorem is thus completed. 

□ 

6.2 The distance problem with mean curvature condition 

We shall now consider the case where A = A(x,k(x)), with k(x) being the 
mean curvature of dfl at the point x. The standard Bernoulli problem 
for harmonic functions (case p = 2) with non-constant gradient condition 
depending on the mean curvature has been studied in |1CH l53j [5l] through 
a variational formulation for the functional 

/ \Vu\ 2 dx + Per({v > 0}). 
Jd 

Here we follow the sub/supersolution approach. 

We first recall the definition of mean curvature in the viscosity sense: 
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Definition 6.2. For U a bounded convex domain, one may define the (in- 
terior) mean curvature of dU in the viscosity sense as follows: assume that 
€ dU and that the interior normal v to d£l at the origin is in the direction 
of the e-axis. We define 

K(dU)(0) := inf k(S a )(0), 
Aeu 

where Sa = {(x,e) : e = (Ax,x)} and il is the set of all symmetric matrices 
A such that the set Sa (the graph of a quadratic polynomial) locally touches 
dU at from inside. 

If the set 11 is empty (i.e., no paraboloid touches dU at from inside), 
we define 

K (dU)(0) = oo. 

We also define the exterior mean curvature 

Definition 6.3. Let U be as in the previous definition: We define the ex- 
terior mean curvature K.(dU)(0) as 

K(dU){0) := su P k(5 j4 )(0), 

where Sa = {(x, e) : e = (Ax, x}} andil is the set of all symmetric matrices 
A such that the set Sa (the graph of a quadratic polynomial) locally touches 
dU at from outside. 

If the graph of any quadratic polynomial locally touches dU at x = 
from outside, we say that 

K(dU)(0) = oo. 

Clearly for any point xq € dU where dU G C 2 in a neighborhood of xq 
it holds that n(dU)(xo) = n(dU)(xo). 

In this section we study the exterior discrete Bernoulli problem (P' E ) with 
distance condition depending on the mean curvature. More precisely, given 
a convex bounded domain K with dK £ C 1,1 we prescribe the condition 

(6.5) dist(s,r,) = ———=:\~(x), x£dn, 

K{dQ)(x) K ' 1 

with I £ (0, 1). We do not include I = 1 since it is not an admissible value 
for the case when K is a ball, as explained in Lemma 16.71 
The main result is summarized in the following theorem. 
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Theorem 6.4. Let K C Mr , N > 2, be a convex open bounded set and 
1,1 € (0, 1). There exist a function u € C 1,a (yi) and a convex open bounded 
domain Cl C Mr 1 , £1 D K, with dQ 6 C ' 1 which solve the following problem: 



(P'k) 



' A p u = in CI \ K, u = 1 in K, u = on <9f2, 
dist(x,Ti) > A k/ -(x) for allx G <9Q, 
dist(x,Ti) = X K j(x) for allx € d£l K , 



where X K j(x) is defined in \6.5\) and d£l K is the set of all points x € <90 for 
which at least one of the following assumptions is satisfied: 

(Al) x is an extremal point of Q and d£l is C 2 -smooth at x, 

(A2) x is not an extremal point, but it can be written as a linear combination 
of extremal points satisfying (Al). 

Moreover the set d£l G C 0,1 does not contain points x of the form 

(A3) x is an extremal point of Q and K(d£l)(x) = k(8£1)(x) = oo, 

(A4) x is not an extremal point, but it can be written as a linear combination 
of extremal points satisfying (A3). 

Recall that extremal points are defined in Definition 12.101 We first show 
a technical lemma that will be used later: 

Lemma 6.5. Let O n be a decreasing sequence of convex bounded domains 
converging in Hausdorff measure to Q. Let xq be a point on d£l; it holds 

(6.6) liminf K(dQ. n )(x n ) < K(dft)(xo), 

{x n } 

where the liminf is taken for all sequences {x n } satisfying x n 6 dVL n and 
x n ->■ x . 

Proof. First note that if K(d£l)(xo) = +oo inequality (|6.6p holds. Let xq be 
a regular point with curvature < k$ := k(8Q)(xq) < +oo. Without loss 
of generality assume that xq = 0. For any e > 0, there exists a paraboloid 
Pko+e with curvature kq + e at the point x = that touches Oil from inside 
at in a neighborhood U e 3 (see Figure [2]). 

Shift now the paraboloid Pk Q + e by a small positive distance 5 in the 
direction of the outward normal u to dVt at 0; denote this shifted ball by 
-ffc +£' Clearly, if 5 is small enough, P% 0+e H dfl C U £ . Since VL n are convex 
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Figure 2: 



sets and {dfl n } converges to 8ft as n /* oo, there exists at least a set 
Q ns , with ra,5 dependent on 5, such that O U (-Pfc 0+e H U e ) C 0„ 5 and <9r2 n<5 D 
dPfc +e = {a^}; moreover P^ 0+e touches tangentially <90 ni5 at the point x na . 
Inequality (16. 6|) follows since e and <5 can be taken arbitrarily small. □ 

The supersolutions are given by those such that 

(6.7) dist(a;,rO> X K j(x), 
and hence we define the set of all supersolutions as 

B := {Q C M. N , convex, bounded, K d Q : uq satisfies (|6.7p |. 

Note that if xo £ <9f2 is a point where k(8CI)(x) = +oo then the distance 
property (|6.7p does not impose any restriction. One may define analogously 
A (and Aq) to be the set of subsolutions (strict subsolutions) through con- 
dition 

(6.8) dist(s,r { )<(<) X K j(x). 
Lemma 6.6. // Qi, Q 2 £ B, then fii n 2 £ B. 

Proof. The proof is a simple consequence of the fact that 

dist(x , {un in n 2 = I}) > max {dist(x, {u^ = I}), dist(x , {un 2 = I})} ■ 

□ 
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Lemma 6.7. The sets B and A are non empty. Moreover, any super solution 
is greater or equal than every admissible subsolution. 

Proof. A possible supersolution for our problem is the function ur con- 
structed in Lemma 13.21 this depends on the fact that there exists R big 
enough such that 

dist(x, Ti) > IR Vx G dB R (0) and I G (0, 1). 

Note that if K is the ball B r (0) and Ct(|x|) the capacitary function in the 
annulus Br(0) \B r (0), condition (|3.ip becomes 

i-ij)'i>i. 

Clearly there exists some finite value for R (dependent on I) that satisfies 
the inequality above if I € (0, 1) but none if I = 1. 

For the subsolution consider the function u £ constructed in Lemma [ 
it holds 

dist(x, {u e = /}) < — , for all x € {u £ = 0}, 
a 

and the desired property (I6.8P is satisfied by choosing 

Q 

e < 



max ie{llfl= i_ E} k(x) ' 

where k(x) is the mean curvature of the level set {ur = 1 — e}. Note that 
max ie{a fl =i-E} K ( x ) < +°° since this level set is a analytic surface (see 
Proposition 12. 8p . 

Lastly we show that any supersolution is larger than every subsolution. 
We proceed as in Lemma 13.31 by assuming that there exists a subsolution, 
denoted by (05,7/5) which is not everywhere larger than a subsolution, 
denoted by (Q s , u s ). We rescale the subsolution until d£l e s touches d£ls at a 
point x G dn £ s n dn s . Then 

K(da £ s )(x ) > K(dn s )(x ). 

Using strong comparison principle we arrive at the following contradiction: 
> dist(xo, {u £ s = I}) > dist(xo, {us = I}) > 



This concludes the proof of the lemma. □ 
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Lemma 6.8. Let {£l n } be a decreasing sequence of bounded convex domains 
in B, and define Vt := Interior(nO,k) ■ Then the p-capacitary potential of 
Q\K , denoted by uq, satisfies (|6.7p and thus, f2 G B. Moreover, there exists 
a positive constant cq such that k(8Q)(x) > cq for all x G Oil. 

Proof. Since each domains f2 n is convex, we can start as in the proof of 
Lemma 13.51 Let u n be the p-capacitary potential of £l n \K. The sequence 
{u n } converges in C 1,a norm to a p-harmonic function u. Convexity imme- 
diately implies that u = outside CI, as it was done in (|3.3p . In addition 
it holds that the set {u n = 1} converges to the set {u = 1} and that d£l n 
converges to dCl in Hausdorff distance. 

We show first that K <s fl: the existence of a subsolution (found in 
Lemma l6,7p implies, using Lemma |2.5( that |Vti n | < M for a constant M 
independent of n. Consequently, for any n G N, it holds 

(6.9) dist(x, {u n = 1})>1-, VxS dtt n . 



This implies that K d Q,. Unfortunately inequality (|6.9p does not give any 
information on uniform bounds from above for n(d£l n ). 

Instead we show that k(8Q) is uniformly bounded from below by a pos- 
itive constant. By construction we have that fi n C do for every n. Then, 
for every point x n 6 dU n , 

(6.10) ' 1j — r < dist(a; n , {u n = I}) < diam(fl ), Mx n £ d£l n . 

K[U\l n )[X n ) 

In particular d$l n cannot have flat parts. Moreover, if we take a sequence 
x n € dVL n such that x n — > xo as n oo with xo G dfl, from (I6.10p and (|6.6p 
(maybe after passing to a subsequence) we get 

< diam(Qo), Vxq G dtt, 



K,(d£l)(xo) 



which shows that d£l cannot have flat parts either. 

Two things are now to be proven: (i) the function u is the p-capacitary 
potential in Q\K and (ii) the function u satisfies (|6.7p . Suppose that Q := 
{supp(u)} C 0; then there exists a point xq G dCl\dQ with a neighborhood 
of radius e, denoted by B e {xq), such that VL U B e (xq) C ft. Let P G {u = 1} 
be the closest point of {u = 1} to xo and cover the segment xq,P with a 
finite number of balls if radius e. Note that this is possible because the 
segment xq,P has a non-negative finite length whose bounds are given in 
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(|6.9p and (|6.10p . By mimicking the argument in Lemma 13.51 one can show 
that Harnack inequality leads to a contradiction. Hence (i) is proven. 

Consider now a point xq £ d£l and any sequence {x n E d£l n } such that 
x n — > xo as n /* oo; it holds 

dist(x n , {u n = Z}) — > dist(xo, {u = I}), as n / oo. 

On the other hand, taking limsup n ^ +00 on both sides of (|6.7p and using 
(|6.6p we obtain 

lim dist(x n , {u n = Z}) > limsup — - > 



K(<9Q n )(x n ) k(<9O)(x )' 

The proof is completed. □ 

Proof of Theorem 16.41 It remains to prove that the p-capacitary 
function (un,£l) constructed in Lemma I6TS1 satisfies the distance property 

(6.11) dist(x,r,) = X K j(x), 

for every x that satisfies either (Al) and (A2) and that <9f2 does not have 
points that satisfy (A3) nor (A4). 

Suppose, by contradiction, that there exists an extremal point xq £ dVl 
such that 

dist(x ,W = Z})> K( J )(xo) + J, 6>0. 

Then there exists a neighborhood U of xo (depending on 5) such that for 
any x 6 U H 9fi it holds 

(6.12) dist(l , { « = 0) >_l_ + |. 

We first assume that dQ is C 2 in a neighborhood of xo, i.e. assumption 
(^41) is satisfied. We apply translations and rotations so that xo = and 
choose coordinates x = (x',x n ) in such a way that Q C {x n > 0} and 
parameterize the surface dQ near xq as a graph x n = (p(x'), x' € M^ -1 , such 
that ip(0) = 0, V</?(0) = 0. Defined a new function <p\(x') := (1— e)i/3(x')+£ 2 . 
For e small enough the intersection set V £ := {(x',x n ), \ ip{x') = ipi(x')} is 
strictly contained in the small neighborhood U. Consider the new function 
(p := max{9?, ip\} and define the set 

n = (n n u n {(a/, x n ) | x n > ^(x')}) u n ^ c )- 
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By construction it holds that O £ fi; moreover the mean curvature of 
can be estimated in terms of the mean curvature or <9f2 at xq, except maybe 
at the points in V £ n dQ, in which it may be infinity. Hence let u be the 
p-capacitary potential of Cl\K. Comparison principle and (16.12|) yield, for e 
small enough, 

15 1 
dist(xi, {u = I}) > dist(a?i, {u n = /}) > + - > , 

k{oU)(xq) 2 k(oU){xi) 

for x\ € 90 n J7. Note that we may obtain the same conclusion of x\ £ V £ 
even if the curvature at those points is ninfinity. The existence of such 
contradicts the minimality of 0. 
Let us assume now that the extremal point xq satisfies (A3). Inequality 
(|6.12p reduces to 

dist(x, {un = I}) > -. 

We perturb Q to in the following way: consider, as before, dSl written 
as the graph x n = ip(x'). By hypothesis, K(dQ)(xo) = oo, so for each 
constant M > 0, there exists a paraboloid m(x') with curvature equal M 
at xq that touches <90 from outside at xo- Define a new function fiix') := 
(1 — e)m(x') + e 2 and consider (p := max{(^,y?i}. Define the set Cl as in 
the previous case; contradiction follows similarly. Hence in such point xq we 
have that dist(xo, {uq, = I}) = 0; but this is not possible because we have 
uniform gradient bounds for uq (see Corollary 12. 5p . 

Next we consider the case of xq not being an extremal point of d£l: 
Lemma 12 in [5l] shows that ^TgQjT^ is a locally concave function if is 
convex. Then we may proceed as in (16.2(1 and following arguments in order 
to conclude that property (|6.1ip is satisfied for each points of the form (A2) 
and that dVt does not contain points of the form (^44). □ 

Remark 6.9. Open problems: 

(i) Full regularity of d£l for (P' K )- For the classical Bernoulli problem full 
regularity of dfl is obtained via the well known regularity theory for 
almost minimal surfaces (see It is an interesting open question 
to show that such theory can be applied also for (P'x)- 

(ii) One may also consider the problem with more general distance func- 
tions \ K f= X(x,k(x)) in (P' K )- 
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6.3 Two- and Multi-phase problems 

The multi-phase Bernoulli problem (c.f. [UEJHS] for p = 2, [B] for general p), 
has many applications in studying interfaces, for instance. In a similar flavor 
we formulate the multi-phase discrete Bernoulli problem. Fix m € N U {0} 
and let K±,K m +2 be two convex bounded domains such that K\ C K m+ 2- 
Moreover consider a sequence of real numbers — l<aj<l,i = 2,...,m + l 
with Oj > etj+i and continuous functions <?j : (K m+ 2\K\) x M + — s> M + , 
i = 2, . . . , m + 1. We seek a sequence {iQ : i = 2, . . . , m + 1} of convex 
domains such that i^i (s (£...(£ K m+ i d -ftr m +2 solve the multi-layer 
distance problem. Let Uj be the p-capacitary potential of the set Ki + \\Ki, 
i.e., solution of 

{A p = in \ Ki, 

Ui = ^ on dKi, 
Ui = a i+1 on 8K i+1 . 

The functions U{ must satisfy the following (nonlinear) joining conditions: 
for < I < min{|aj — aj+i|} the level sets 

r* +z := {ui = a i+1 + 1}, rj] 1 := {u i+1 = a i+1 - 1} 

satisfy 

(6.13) Fi (x,dist(x,r^. 1 ) J dist(x,r^ 1 )) = 0, for all x 6 8K i+1 , 

for given functions Fi : K x M + x M + — > M, i = 1, . . . , m. 

We assume that Fi(x,p,q) is a continuous function, and that Fi(x,p,q) 
is strictly decreasing as a function of the variable p for all g. Then via the 
implicit function theorem one can write the joining condition (16. 13f) as 

dist(x, r\.,) = g (x,dist(x,T i + 1 )) , on dK i+1 , 

for some functions g : K x M + — > M + . 

We assume the following initial hypothesis on the functions g%, i = 
l,...,m: 

(HI) gi is continuous and bounded from below by a positive constant. 

(H2) gi is non-decreasing with respect to the second argument. 

(H3) The function x i— > gi (x, q(x)) is a concave function whenever q(x) is a 
concave function. 
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(H4) For any given value yo > 0, % = 1, . . . ,m, there exist constants < 
ci < c 2 such that c\ < gi(x,y)/y < c 2 , uniformly for all x £ K and all 
2/ > Vo- 

Although condition (H4) may seem artificial at first, it is the classical 
hypothesis that assures convexity for the level sets of the solution to the 
Bernoulli problem (c.f. [IS]). The classical example is 

9{x ' q)= (0(^ + ^)1/"' 

for a function a such that l/a(x) is concave. Existence, uniqueness and 
convexity of solutions for a general g is an open problem. 

The main results of this sections are in the following two theorems: 

Theorem 6.10. Let K U K 3 C R N be two convex bounded domains such 
that K\ <<= K 3 . Then there exists a convex domain K 2 such that 

K x t=K 2 m K 3 , dK 2 G C 1 ' 1 , 

and the p-capacitary potentials u\ and u 2 of the sets K 2 \K\ and K 3 \K 2 , 
respectively, i.e., solutions of 

{A p ui = in K 2 \Ki, ( A p u 2 = in K 3 \K 2 , 
m = 1 on dK u I u 2 = on dK 2 , 

u\ = on dK 2 , [ u 2 = — 1 on dK 3 , 

satisfy the nonlinear joining condition 

(6.15) dist(x, {u\ = /}) = g (x, dist(x, {u 2 = —I})) on dK 2 . 
Theorem 6.11. The multi-phase version of Theorem \6.10\ holds. 

The proof of Theorem 16.101 already contains the fundamental ideas for 
m > 2: indeed Theorem 16.111 follows easily from Theorem 16.101 and the 
uniform separation results of Theorem 16.151 

Consider the class of convex sets 

C := {SI C K , convex bounded domain : K\ iOg i^} , 

and u\, u 2 the p-capacitary potentials of K\K\ and K 3 \K, respectively. 
Define the set of subsolutions and supersolutions as 

A := <E C I dist(x, {ui = I}) < g (x, {u 2 = -I})} , 

(6.16) B:={n€C\ dist(x, {ui = I}) > g (x, {u 2 = -/})} , 



40 



for x € dK2- 

We next show that the set B is closed and stable under intersections: 
Lemma 6.12. 7/0,Oe B, then nnCleB. 

Proof. Let m, U2, ui, U2 be the p-capacitary potentials of Q\Ki, K^\Q, (l\K\ 
and K^\^l respectively. By construction n is a convex domain. Next, 
let v\ and v-i be the p-capacitary potentials of (J) n f2)\i£i and i^3\(f2 PI Ci) 
respectively. Let x € c?(£7 n O); without loss of generality we assume that 

x G on. 

By comparison principle, v\ < min{ui,ui}. It holds 

(6.17) dist(x, {v\ = I}) > dist(cc, {u\ = I}) > g (x, dist(x, {u2 = —I})) , 
since f2 € B. 

On the other hand, by comparison principle we have that V2 > max{u2, U2}; 
thus 

dist(x, {V2 = —I}) < dist(x, {112 = —I})- 
Since g is non-decreasing with respect to the second argument, it holds 

(6.18) g (x, dist(x, {V2 = —I})) < g (x, dist(x, {112 = —I})) ■ 
Formula (I6.17P together with (I6.18P gives that 

dist(x, {v\ = I}) > g (x, dist(x, {V2 = —I})) , 
which implies that finfi £ B. The proof of the lemma is completed. □ 

Lemma 6.13. Let Hi D O2 D ■ ■ ■ be a decreasing sequence of convex do- 
mains in B. Then := Interior^Q,^) 0, belongs to the class B. 

Proof. Let u\ and u\ be the p-capacitary potentials of £l\\K\ and K^Q^, 
respectively. By standard arguments one can show that u\ and u\ converge 
in C a to vi and V2, the p-capacitary potentials of Q\K\, Ks\£l. Moreover 
the convergence is C 1,a away from the boundary sets. 

To conclude the proof of the lemma, we need to show that dist(x, {v% = 
I}) > g (x, {v2 = —I}) for x € d£l: the proof mimics the one of the one-phase 
case (see Lemma 1331 □ 

Lemma 6.14. The classes A, B are nonempty. 
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Proof. Consider the solution u of the p-capacitary problem in K^\Ki\ 

A p u = in K 3 \K[, 
u = 1 on 8K\, 
u = — 1 on 8K 3 . 

For any e € (—1, 1), set K e := {u = e}, and define the function 



(6.19) 



u £ (x) :-- 



. . u(x) OL „ r . . T . 

1 — a 

, . u(x) — a , \ 

U 2 ,e( x ) = "~ T~ fOT X G ^3\^ £ - 

1 + a 



Then one may check, using assumption (H4) that K £ € B for e close to — 1, 
and that K £ € A for e close to +1. □ 

Proof of Theorem 16. 101 Let K2 be the minimal set in B; K2 is convex 
and well defined by the previous lemmas. Let u±, 112 be the p-capacitary 
potentials of K2XK1 and K 3 \K 2 , respectively. Then we just need to check 
the joining condition (|6.15p . Suppose by contradiction that the property 
(I6.15P fails. Then there exists xq € dK 2 such that 

dist(xo, {u\ = I}) > g (xq, dist(x, {u 2 = —I})) + £■ 

By continuity, the same inequality holds (maybe with e/2) in a whole neigh- 
borhood of xo- If xq is an extremal point of K2 the same argument as in 
Section 16.11 gives a contradiction. 

However, if xq is not an extremal point but is a finite linear combination 
of extremal points, we need the following additional argument: for x £ 8X2, 
x 1 y dist(x, {iti = I}) is a convex function and x h-> dist(x, {112 = —I}) is 
concave. By the concavity assumption (H3) of g, we know that the function 

G(x) := dist(x, {tti = I}) — g (x, dist(x, {11,2 = —I})) 

is also convex. Since K2 € B, we know that G > 0. But we have just 
shown that G vanishes for all extremal points of K2. Then, if xo is a finite 
linear combination of extremal points, then we must have G(xq) = 0, which 
concludes the argument and also the proof of Theorem 16.101 

For the regularity of the free boundary, just note that if g is bounded 
from below by a positive constant, then 8X2 satisfies the interior ball con- 
dition. Since K2 is convex, we may conclude that 8K 2 G C 1 ' 1 . □ 

We finish this section with a uniform separation result, in the spirit of 
Theorem 5.1 from [6] for the classical Bernoulli problem: 



42 



Theorem 6.15. Let H be the set of configurations (K\, K 2 , K%) such that 
K\ , K 2 , K% are convex bounded domains satisfying 



• Ki (s K 2 K 3 , 

• K\ satisfies the interior ball condition for radius > t\, 

• K 3 satisfies the interior ball condition for radius > r 3 , 

• K 2 belongs to the set of supers olutions B from (|6.16[) . 
Then there exists a value n = n(R,ri,r 3 ) such that 

dist(dK 1} dK 2 ) > •qdist{dK 1 ,dK 3 ), 
uniformly for all (K\, K 2 , K$) € H. 

Proof. Without loss of generality, assume that dist(dKi,dK 3 ) = 1, other- 
wise the result follows by rescaling. Let u be the solution of the p-capacitary 
problem (|6.19p and set 

a := sup{u(x) : x 6 dK 2 } £ (—1, 1). 

We claim that there exists a uniform ao such that for every admissible 
configuration we must have a < a^. For the proof of the claim, it is enough 
to restrict to configurations that have a € (0, 1). Let u\,u 2 be as in (|6,14p . 
Consider the functions 

, u(x) — a . . u(x) — a 

ui,ai x ) '■= — i 5 u 2 , a {x) := — — 

1 — a 1 + a 



in the set K 3 \K\. Then u± = 1 = u\ a on dK\, while ui )Q < = u\ on dK 2 . 
Comparison principle gives that u\ a < u\ in the set K 2 \K\. Similarly, since 
u 2 = — 1 = u 2ta on dK 3 and u 2)a < = u 2 on dK 2 , we arrive at u 2<a < u 2 
in the set K 3 \K 2 . 

Choose xq £ dK 2 such that u(xo) = a. Because K 2 is a supersolution, 
we must have 

dist(x , {tii = 0) > 9(xo, dist(x , {u 2 = -/})). 
On the other hand, since we had that u\ a < u\ in K 2 \K±, 

dist(xo, {tii = 1} < dist(xo, {iii, Q = /}) = dist(xo, {ii = /(l — a) + a}). 
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In addition, the relation U2 a — u 2 in K%\K2 gives 

dist(xo, {u2 = —1} > dist(xo, {^2,0 = — I}) = dist(xo, {u = —1(1 + a) + a}). 

By hypothesis, g is a non-decreasing function, so the previous three inequal- 
ities give the relation 

(6.20) dist(xo, {u = 1(1 — a) + a}) > g(xo, dist(xo, {u = —1(1 + a) + a})). 
Taking into account that u(xq) = a, we may now estimate 

1(1 — a) 1(1 — a) 
dist(aj , {« = K 1 - «) + «}) < ■ riw ^ 



inf|Vn| mi 

for some mi depending on ri, where we are using Lemma 12.61 to bound the 
gradient from below. Moreover, 

dist(x , {u = -Z(l + q) + a}) > ^ — y , > -i 

sup|Vu| m,2 

where m2 = ^(t^) is given in Corollary 12.51 Formula (|6.20p and the mono- 
tonicity of g yield to 

1(1 -a) ^ ( 1(1 + a) 
>g[xo, 



mi \ m2 

By hypothesis (HA) on g it holds 

(0.21) <<IZ£>> ei .'< 1 + <»> 



mi 771,2 

for ci depending on 7712- 

Note that mi and 7712 depend on the distance dist(9-Ki, dK%) = 1, and 
on the initial constants ri,r%. But then, (|6.21|) already implies that a < 
for some «o- The claim is proved. 

Finally, to complete the proof of the theorem, one uses Lemma 5.3 in 

m- □ 

Remark 6.16. If we relax the hypothesis on the domains Ki,K% to only 
a interior cone condition, then we can use the uniform C@ estimate from 
Corollary \2.4\ instead of the gradient estimate, and still obtain a uniform 
separation result. 
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6.4 A Brunn-Minkowski inequality 

Given two domains fii and U2 C K^, we define their Minkowski linear 
combination fit as 

n t := (1 -t)fi + te[o,i]. 

Notice that if fio, ^1 are convex sets, so is fit. 

An upper semicontinuous function u : 4lU {±00} is said to be 
quasi-concave if it has convex superlevel sets, or, equivalently, if 

u((l — t)xQ + tx\) > mm{u(xo), u(x\)}, for all t G [0, 1], xq,x\ G M. N . 

If it is defined only in a proper subset fi in M^, we extend u as —00 in 
IR^fi and we say that u is quasi-concave in fi if such an extension is quasi- 
concave in Mr. In an analogous way, u is quasi-convex if — u is quasi-concave. 
Obviously, if u is concave (convex), then it is quasi-concave (quasi-convex). 

Consider uo,u± two upper semicontinuous functions defined in fio,fii C 
K , respectively, and let t G [0,1]; the Minkowski linear combination of 
uq and U\ is the upper semicontinuous function u\ whose super-level sets 
C? := {tij > 1} are the Minkowski linear combination of the super- level 
sets Cf\ dp of uq,u\, respectively, i.e., 

and 

u*(x) = sup{/ : x G df^}. 

The notion of quasi-concavity has already been used in the study of the 
classical Bernoulli problem in the papers [151149]. for instance. Here we plan 
to extend those results to the distance problems (Pe) an d (Pi)- It turns 
out that the proofs seem to be very well adapted for the case of distance 
between level sets. 

Our first proposition deals with the exterior case: 

Proposition 6.17. Fix I G (0,1). Let Kq, K\ be two bounded convex do- 
mains in M , and Xq, X± two given positive constants. For t G [0, 1], define 
their Minkowski sum 

K t := (1 — t)K + tKi, and X t := (1 - t)X + t\ x . 

Denote by (uo,Qq), («i,fii), (ut, fit) the solutions of the exterior distance 
problem (Pe) for given data (Kq,Xo,1), (K±,Xi,l) and (K t ,Xt,l), respec- 
tively. Then 

(6.22) (1 — t)fio + ifii D fit- 
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Proof. Let := (1 — t)£lo + tCl%; Consider the Minkowski sum u\ of the 
functions uq and u\. Then for every level set s G [0, 1], we have that 

{u* t >s} = {\- t){u >s} + t{ui > s}. 

Let x G d£lt and y G {u$ = I}. Then there exist xq G 8VLq and x\ G d£l\ 
such that x = (1 — t)xo + ixi and as a consequence, 

dist(x,y) > (1— t)dist(xo ; {^o = l})+tdist(xi, {u\ = /}) = (1— t)\Q+t\\ = \ t . 

This implies that Vtt G B\ t . 

Since the solution of the problem (Pe) for given initial data (K t ,X t ) is 
given as the minimal set in the class B\ t , we arrive at Qt D fit, as desired. □ 

Remark 6.18. We conjecture that equality in (|6.22p holds if and only if Kq 
and K\ are homothetic. 

We consider now the interior counterpart: let I G (0, 1) and be £1 a 
convex domain in H N . Denote by A(J1) := Xn j7n ax the Bernoulli constant 
for the (interior) distance problem defined in (|4.4p . 

Proposition 6.19. Fix I G (0,1). Let Qq, Q\ be two bounded convex do- 
mains in H N , and Ao < A(f] ), Ai < A(fii) two given positive constants. 
For t G [0, 1], define 

n t := (1 -t)U + and X t := (1 - t)A + iAi. 

Denote by (uq,Kq), (ui,K\), (ut,K t ) the solutions of the interior distance 
problem (Pi) for given data (£lo, \o,l), (Qi,Ai,Z) and (Tit, At, Z), respectively. 
Then 

(6.23) (1 - i)lfo + t-fiCi C i^t- 

Proof. By definition of the Bernoulli constant (|4.4p , there exist convex sets 
Kq <<= SIq, K\ (g fli, whose p-capacitary potentials iio and ui, respectively, 
satisfy dist(xj, {itj = /}) > Aj, for all X{ G dKi, i = 0, 1. 

Let uj? be the Minkowski addition of the quasi-convex functions uq and 
iti; u* is a quasi-convex function whose sublevel sets are the Minkowski 
linear combination of the corresponding sublevel sets of uq and u\, i.e. 

{ u * < s } = (1 - t){n < s} + t{m < s} for all s G [0, 1]. 
We define K t := {u* t = 0}; note that K t = (1 - t)K + tK x . 
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Then by a purely geometrical argument we can show that if x € dKt 
and y S {u* = I}, then we may find xq G 9Kq, x\ € such that 

dist(x, y) > (1 — t)dist(xo, {no = /}) + tdist(a;i, {ui = /}) > (1 — i)Ao + iAi. 

Hence for all x 6 dA't it holds 

dist(x,{iij = /}) > Af. 

The inequality above implies that the set VL t belongs to the class B\ t . Since 
the set Kt, solution of the interior distance problem (Pi), is constructed as 
the largest set in the class B\ t , then one automatically obtains that 

{l-t)Ko + tK 1 = K t C K t . 

□ 

As a consequence, we obtain a Brunn-Minkowski inequality for A(f2) in 
the spirit of the (Newtonian) capacitary inequalities of [16] : 

Corollary 6.20. Fix I € (0, 1). Let f2o>^i be two bounded convex domains 
in R^. For each t € (0, 1), define 

n t ■= (i-t)n + tni. 

Then 

(6.24) A(fit) > (1 - t)A(fi ) + tA(Ox). 

Conjecture 6.21. We conjecture that equality in (|6,24p holds if and only 
if S7o and £l\ are homothetic. One may also compare the Bernoulli constant 
A(O) to the Bernoulli constant of a ball with the same mean width, and 
obtain (sharp) isoperimetric inequalities. Another open problem is whether 
uniqueness of solution for A = A(il) holds. 

7 Concluding Remarks 

7.1 An alternative approach for Numerical computation 

One of the major problems in designing numerical algorithms for the (inte- 
rior and exterior) classical Bernoulli problem, as well as for the formulation 
presented in this manuscript, is related to the fact that the set f2 (for the 
exterior case), or K (for the interior case) is unknown and chosen by the 
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solution. Consequently, to impose numerically the gradient or the distance 
condition to the solution itself at the unknown boundary becomes a non- 
trivial task. 

Consider for a moment the exterior discrete Bernoulli problem (the in- 
terior can be treated similarly) 



{A p u = in Q, \ K, u = 1 in K, u = on dQ, 
dist(x, {u = \oj}) = A for all x E dQ, 

and the classical one 

A p u = in £1 \ K, u = 1 in K , u = on dtt, 
iViil = uj for all x € d£l. 



(Pe) 



(Pb) 



The numerical implementation of a distance condition between two level 
sets of a p-harmonic function might be easier than the implementation of 
a gradient condition on a level set. Moreover we have rigorously shown in 
the previous sections that solutions of (-Pe) converge to solutions of (Pb) as 
A ->■ 0. 

Based on these two observations, we propose the following algorithm 
to numerically approximate solutions to problem (-Pe) an d consequently 
solutions to (Pb)- Fix the set K C M> N and (u n _i, f2 n _i) solution of (Pn) 
with f2 = Q n -i: 

1. Let P>\(x) be the ball of radius A and centrum in x; define 

n n := {n n „! > \cj} U J S aO) 

with r n _i = {x £ O n _i, I u n -\(x) = Xco}. 

2. If fi n _i = Q n the couple u n -i) is a solution for (Pe)- If 
f2 n _i ^ fi n solve (P^) with f2 = O n . 

3. Let (u n ,£l n ) be a solution of (Pn) with 17 = £l n . Restart from step 1. 
for n n + 1. 

Hence a numerical algorithm for (Pb) consists on solving (Pe) with the 
scheme {1., 2., 3.} for a decreasing sequence of values for A. 

We remark the reader that estimates showing convergence and stability 
of algorithm {l.,2.,3.} are not provided in this manuscript but are very 
interesting open questions. 
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7.2 Open problems 

Let K m £1 be two bounded domains in R N . Let u G C 2 (Q\~K) n C(H) be a 
classical solution of the Dirichlet problem 



where F(x,t,p, A) is a proper and (degenerate) elliptic operator defined 
on M. N xMx M. N x S N (here 5^ denotes the set of real symmetric n x n 
matrices). If both K and f2 are convex, the natural question to ask is whether 
all the level sets of u are convex, i.e. whether u is a quasi-convex function. 
Without suitable assumptions on F, the answer can be negative (see [55]) for 
instance. Before dealing with the classical Bernoulli or the discrete version of 
it presented in this manuscript, one should aim to find sufficient conditions 
on F which guarantee that u is a quasi-convex function. An answer to 
this question for the p-Laplacian operator can be found in Lemma 12.81 and 
Lemma 15.31 for respectively the convex and star-shaped domains. Note that 
for a general F, the problem was considered in [58]. 

Another interesting open question is whether any solution to (Pe) (° r 
to (Pi)) satisfies the so called normal vector property: a solution (u,Cl) to 
(Pe) satisfies the normal vector property if given any point xq 6 Oil the line 
perpendicular to the boundary dQ at xq intersects the convex hull of K, as 
in Figure ??. Similar questions have been addressed for reaction-diffusion 
equations in the context of front-propagation (Jones' Lemma, see |39|. I12j). 

Finally, one may also consider cases in which the condition (jl.ip is re- 
placed by 



or, even more general, by the distance between any two levels sets. Such 
problems fall outside the scope of this paper and are left as open questions. 

A still open and tantalizing problem concerns the formulation of the dis- 
crete Bernoulli problem considered here via a suitable variational functional. 
The next question would then be to formulate suitable isoperimetric prob- 
lems, in the spirit of |24j . 
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